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Abstract 

For the stability of agglomerates of micron sized particles it is of considerable 
importance to study the effects of tangential forces on the contact of two particles. 
If the particles can slide or roll easily over each other, fractal structures of these 
agglomerates will not be stable. We use the description of contact forces by Johnson, 
Kendall and Roberts, along with arguments based on the atomic structure of the 
surfaces in contact, in order to calculate the resistance to rolling in such a contact. 
It is shown that the contact reacts elastically to torque forces up to a critical bending 
angle. Beyond that, irreversible rolling occurs. In the elastic regime, the moment 
opposing the attempt to roll is proportional to the bending angle and to the pull-off 
force P c . Young’s modulus of the involved materials has hardly any influence on 
the results. We show that agglomerates of sub-micron sized particles will in general 
be quite rigid and even long chains of particles cannot be bent easily. For very small 
particles, the contact will rather break than allow for rolling. We further discuss 
dynamic properties such as the possibility of vibrations in this degree of freedom 
and the typical amount of rolling during a collision of two particles. 


§ 1 . Introduction 

The coagulation of small solid particles is a process of great importance in a wide 
range of fields such as material sciences (sintering) atmospheric sciences (aerosols) and 
also astrophysics (the formation of planets; see, for example, Weidenschilling, Donn 
and Meakin (1988)). Many studies of coagulation have been carried out both on the 
experimental side and on the theoretical side. Coagulated aggregates have structures 
that vary from quite compact to open fractal. Numerical simulations and also theoretical 
arguments (see, for example, Ball and Witten (1984)) have shown that the density of 
aggregates reflects the formation process, in particular the type of trajectory of the 
approaching particles (ballistic or random walk). A further important factor is the 
difference between particle-cluster aggregation and cluster-cluster aggregation. In 
particular, the addition of single particles to the cluster gives rise to quite compact 
structures as compared to cluster-cluster aggregation. 

However, examination of pictures of (especially low density) aggregates shows that 
the stability of the adhesive joints is also an important factor for the structure (Kantor 
and Witten 1984). Often, the aggregates consist mainly of quasi-Iinear chains of 
particles with no more than two contacts per particle. Only strong joints can survive 
collisions with further particles without tearing, and only if the joints can resist 
tangential forces will the chains continue to be linear. If, for example, rolling of one 
grain.over the other were easy, the chains would quickly fold and the aggregates would 
become considerably more compact. 
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As far as the breaking of the aggregate is concerned, certainly the most important 
parameter is the pull-off force 


P c = 3nyR , (1) 

which is the force required to separate two spherical particles in an adhesive contact. 
Here 2y is the specific energy of adhesion for the two surfaces and \/R = \/R l + \/R 2 
is the reduced radius of the particles involved. The factor 3 is the value obtained from 
the JKRS theory (see below). For very hard spheres, this factor should rather be 4 
(Muller, Yushenko and Derjaguin 1982). 

When the forces involved are not sufficiently large to destroy an aggregate, it is still 
possible to deform it by bending the particle chains, rolling the spheres over each other 
or sliding the surfaces along each other. The latter is not likely to be too important after 
the formation of the aggregate, since the forces needed to break a contact by horizontal 
forces are usually larger than the pull-off force (Savkoor 1992). However, what should 
be much easier is bending these chains, and eventually moving the contacts by rolling 
the particles over each other. 

To get a handle on this, it is important to know which forces oppose the attempt 
to roll one body over the other in an adhesive contact. In this paper, we will develop 
a model which invokes both the properties of an adhesive contact between two elastic 
spheres and, to some extent, the molecular structure of the surfaces in order to calculate 
these forces and to estimate the energy dissipated in the rolling process. 


§ 2. The JKRS solution for the contact problem between two 

ADHESIVE SPHERES 

Our considerations are based to a large extent on the solution of the adhesive contact 
problem which was developed by Johnson, Kendal and Roberts (1971) (JKR) and also 
by Sperling (1964). The JKRS solution is based on fracture-like energy considerations. 
The main assumption leading to this solution is that the forces between the two spheres 
are pure contact forces. That is, two grains feel only an attractive forces across the 
contact area, while the surface areas outside the actual contact area are force-free. 
Certainly, such an approximation has its limitation. It was shown by Muller, Yushenko 
and Derjaguin (1979) that JKRS is more appropriate in the case of large soft spheres. 
For small and very hard spheres, the description by Derjaguin, Muller and Toporov 
(1975) is more accurate. We have decided to adopt JKRS for the time being, since the 
assumption of contact forces considerably simplifies the discussion of rolling 
resistance. 

JKR showed that two elastic adhesive spheres make contact over a finite circular 
region, as indicated in fig. 1 . In the case of no external applied forces, the radius of this 
contact area is 


a 0 = 



( 2 ) 


where y is the surface energy of each surface. For different surfaces, y = y , + y 2 - 2y n , 
where y n is the interface energy. R is the reduced radius of the two spheres 
(R 1 = R\ 1 + Rz ') and (£*) 1 = (1 — v\)IE\ + (1 — v\)/E 2 , where E, and v, are 
Young’s modulus and Poisson’s ratio, respectively, of grain The elastic displacement 
of the two spheres along the line connecting the centres of the two spheres (again with 
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Fig. 1 


z 



coordinate origin 

Sketch of the geometry of the contact area between two elastic spheres. The two bodies make 
contact over a circular region with radius a . Inside the contact region, the surfaces interact 
by a mutual pressure p (which in adhesive contact can also be negative). When tangential 
forces are applied, they are transmitted by tangential tractions q in the contact area. The 
contact area is always in the x-y plane. 


zero external force) is given by 



(3) 


where — P c = — 3nyR is the maximum external force that can be applied without 
breaking the contact (the so-called pull-off force). In the following we will sometimes 
use dimensionless quantities, denoted by a hat over a symbol: 


P = 


P 

K' 




(4) 


Here P is the applied external load. When this load differs from zero, both the contact 
radius and the displacement change. The following relations between load, displace- 
ment and contact radius hold for the JKRS solution: 


P = 4 (d 3 - d 3/2 X (5) 

3 = 3d 2 - 2 d m . (6) 


The distribution of pressure in the contact area is given by a superposition of Hertz’s 
and Boussineq’s contact problems: 

p(r, a) = pu( 1 “ f/a 2 ) + p B ( 1 “ fta 2 ) ~ m . (7) 

In JKRS, the two constants are given by 


P h 



P b = 



( 8 ) 


(9) 
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§ 3. Resistance to rolling in a JKRS contact 
Resistance to rolling is a force that opposes the attempt to roll a body on the surface 
of another body. To measure this force we might apply a force Q x at the centre of a 
sphere with radius R\ in the x direction parallel to surface on which the sphere rolls. 

In order to study this effect independently of the sliding process, we assume that 
sliding friction is high enough to prevent sliding entirely. Then, the tangential force Q x 
is balanced by an equal force in the contact area. However, since the two forces are 
applied at different points on the sphere, there is a resulting torque about the y axis 
which is 


M y = Q X R . (10) 

In order to prevent the sphere from starting to roll, this torque has to be balanced by 
an appropriate pressure distribution in the contact area. In order to produce a resultant 
moment (the so-called rolling moment), the distribution would have to be asymmetric 
with respect to the y axis. However, the pressure distribution in JKRS is always 
cylindrically symmetric with respect to the contact point. Since this implies 
p(x,y ) = p( — x ,y), there is no moment associated with this distribution: 

M y = J Jxp(x , y) dx dy = 0. (11) 

Different mechanisms have been discussed in order to explain the observed rolling 
friction at least in the adhesion-free problem. The different possibilities are micro-slip 
at the interface, ineleastic deformation of the involved materials and surface 
irregularities; these have been summarized by Johnson (1989). However, none of these 
mechanisms should be of great importance in the contact of (sub-) micron-sized 
particles, of which adhesive aggregates usually consist. The displacements at the 
surfaces of the particle are usually not much larger than one or a few diameters of the 
atoms in the solid. Therefore, the movement of dislocations, which for macroscopic 
crystalline bodies is the source of inelastic behaviour, will play no role at very small 
particle sizes. Likewise, surface roughness is unimportant on the 50 A size scale of the 
contact area of sub-micron-sized particles. 

However, for very small spheres, the surfaces can no longer be treated as perfectly 
smooth at all scales, since they are made of atoms. Then it becomes clear that the contact 
area cannot be perfectly symmetric. During rolling, new contact at the leading edge of 
the contact will not be made continuously, but rather in steps after a rolling distance 
of approximately one grid constant of the material. Likewise, at the trailing edge, the 
contact is lost in steps. Therefore the pressure distribution will not be exactly 
symmetrical all the time. 

Still, if new contacts at the leading edge are made very smoothly and if the processes 
at the leading and trailing edges are exactly time reversed to each other, no energy will 
be dissipated and the pressure distribution will still be symmetric on average. This 
means that at every instant there may be a non-vanishing torque, but this torque may 
change sign and the average will still be zero. However, it is known from experiments 
with the atomic force microscope (Landman, Luedtke, Burnham and Colton 1990) that 
the approach of a tip to an adhesive surface and subsequent retraction involves a 
hysteresis. Upon close approach, the tip will form an adhesive contact with the substrate. 
Upon retraction, however, the tip will pull out a neck of surface material and contact 
will not be broken until a critical displacement is reached. A similar process occurs for 
a single atom attracted by van der Waals forces to the substrate and held back by 
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Fig. 2 
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Sketch of the mechanism of rolling friction. Top: side-view on an elastic sphere in adhesive 
contact with a rigid plane. Centre: top-view on the contact area. The dotted circle 
indicates the position of the contact point (the projection of the sphere centre). From left 
to right, the diagram shows (a) an initial torque-free state with contact area centred 
around the contact point, ( b ) a state with resulting torque where the contact area is not 
centred and (c) the final and again torque-free state after the contact area has readjusted. 
Bottom: decomposition of the asymmetric contact area into two half-circles of different 
radii. 


spring-like binding forces in the approaching tip (see, for example, the discussion in 
Israelachvili ( 1 992)). Upon approach, the atom will jump into contact with the substrate. 
The energy difference between the two states before and after the jump is dissipated. 
The same process occurs upon retraction between a different pair of states. It has been 
suggested by Tom&nek (1993) that this hysteresis might be the source of atomic scale 
rolling friction. In this paper, we follow up this idea and try to calculate the forces 
involved. 

A possible approach to calculate the rolling friction in an adhesive contact would 
be to look directly at the energy dissipated in the making and breaking of contacts. This, 
however, requires an accurate knowledge of the forces between the atoms and an 
accurate description of the transition between the contact area itself and the surrounding 
parts of the spheres. 

We therefore follow a different approach. The retarded making and breaking of 
contacts means that the contact area will not always be centred around the line 
connecting the centres of the two spheres, but rather will be shifted by a certain amount 
to one side. This shifted contact area leads to an asymmetric pressure distribution, as 
we will show below. 

For now we restrict out terminology to the case of a rolling elastic sphere on a rigid 
plane. However, the case of two elastic spheres can be treated by simply replacing the 
sphere radius R and the elastic constant £* with reduced quantities. 

To describe the rolling process we visualize it as shown in fig. 2. 

(1) We start from a situation where the contact is a JKRS contact centred around 
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the contact point which is the point vertically below the centre of the sphere. 
With two spheres it would be the point on the line connecting the two centres. 

(2) Then the sphere is rolled a small distance over the surface. If d xy denotes the 
angular motion, the contact point moves a distance £ = RS ay . However, we 
assume that the contact area remains behind. Thus the contact area is no longer 
centred around the actual contact point, but is shifted by a distance In this 
state there will be a non- vanishing torque opposing the rolling direction. In § 3. 1 
we will calculate this torque. 

(3) At a certain critical rolling distance £ yield (which should be typically of the order 
of one atom diameter in the material) the contact area readjusts — that is, contact 
is lost at the trailing edge and new contact is formed at the leading edge — so 
that again a symmetrical (torque-free) distribution is formed. 

In the following we will calculate the moment that is produced by the asymmetric 
pressure distribution. 

3.1. The moment associated with a shifted contact area 

To our knowledge, there are no analytical solutions for the problem of a shifted 
contact area. To calculate the pressure distribution and the associated moment we 
proceed as follows: 

(1 ) The true shifted contact area is approximated by a decomposition into two half 
circles of different radii a + <!;, respectively (see fig. 2). 

(2) The half circle x < 0 is one half of a symmetrical contact with contact radius 
a + f and pressure distribution p{r , a + £ , 5 Z ) (see below for the definition of this 
pressure distribution). 

(3) the half circle x > 0 has a smaller contact radius a — £ and a corresponding 
pressure distribution p(r,a — £, d z ). 

(4) The total momentum due to both half circles is calculated. 

The resulting distribution is discontinuous at x = 0. As indicated in fig. 2, the true 
contact area is closely fitted far from the y axis where the contributions to the moment 
are most important. In Appendix 1 we discuss the two-dimensional (2D) case where 
the exact pressure distribution can be obtained analytically. That analysis shows that 
the above approximation gives reasonable results. 

We start from the general pressure distribution for the contact of two elastic spheres 
with a contact radius a . This pressure distribution is the superposition of the solutions 
of the Hertz problem and of the Boussineq problem. The subscripts H and B indicate 
that in the following: 

p(r,a,5 z )=p H (l - ^/a 2 ) 1 ' 2 + . 0 2 ) 

The pressure distribution corresponds to a normal displacement of the spheres 

7 ia 

3z = ^^.(Ph + 2p B ). (13) 

To calculate the moments of one half of the contact area for different sizes for the 
contact area we cannot directly use the JKRS equations since the different contact radii 
a always correspond to different values for S z . For the intended calculation, however, 
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we need contact areas of different sizes but with the same normal displacement 6 Z . We 
have thus to start from eqn. (12). 

For a contact area of a given size, the constant p H has the fixed value = 2aE*lnR , 
since it is this part of the solution that deforms the spherical surface to a flat area. 
Different values of S z can be achieved by changing only the Boussineq part of the 
solution. Therefore we will now replace the Boussineq constant p B with eqn. (13). Then, 
we may keep 5 Z constant and change the size of the contact area. 

An expression for p B can be obtained from eqn. (13): 

Pb(o) = -S z - ip H (a). ( 1 4) 

na 


We insert this into eqn. (12) and get 


p(r, a y d z ) =p H (l ~ P/a 2 ) m + 


(E*fna)d z — pw/2 
(1 — Pfa 2 ) m 


(15) 


As discussed above, we compose the asymmetric contact region of two halves. The 
distribution inside each half is taken to be that of a symmetric pressure distribution with 
contact radius a + c and a — respectively: 


Pain a, 0 = 


j>(r,a + £,S Z ), x<0, 

l p(r y a — £,d z ), Jt>0. 

The moment associated with this distribution is 

r3ll/2 ra + £ rll!2 ra — £ 


(16) 


rjll/2 + £ rTlll ra — C 

My = I I xp(r y a + S z )rdrd(p + I I xp(r, a — S z )rdrd(p. (17) 

JnI2 JO J - Ilf. 2 JO 

With x = rcos q> we may carry out the <p integration and find 

r a + £ r<* ~ G 

M y = — 2 Pp(r , a + £, S z ) dr 4- 2 Pp(r , a — £, <5 Z ) dr 

Jo jo 


= - MX 


( 1 / 2 ) 




* >\ a + £ iy * y, a — £• 

Thus we have to solve integrals of the kind 


(18) 

(19) 


A^!5U = 2 f Mr, sf, 6d d r. (20) 

•>0 

The required integrals are given by eqns. (87) and (88) in Appendix B. The result is 

M < y ][ tf'$ z = 2^/? H “^ 3 + ( 21 ) 


= ±E*st 2 5 z -^p H rf\ (22) 

o 

We insert the correct value of pn (eqn. (8)): 

M[ U2 J,S. = ^ <5 M 2 - \Pc~T ■ (23) 

Z Oq 

The moment of the whole asymmetric contact can now be calculated from eqn. (19). 
We keep S z constant and neglect all terms of higher than first order in £fa. This is 
equivalent to the substitutions srf 4 — » — 8 a 3 £ and stf 2 — > — 4a£: 
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M y » 2{( - E*aS z + 3 P c d 3 ). (24) 

Now we may insert the correct value S z which is the value from Johnson’s solution 
corresponding to a contact area radius a as given in eqn. (6). The result is 

M y = 4P c a 3r2 Z. (25) 

The moment associated with a contact area shifted by a distance £ is therefore 
proportional to the pull-off force and to the shift itself. The order of magnitude of this 
result is what one would have expected, since the forces involved are of the order of 
P c and the force is applied over a level of length £. The dependence on the normal load 
P enters through the factor d 3/2 , which should be close to one in most cases (see below). 

In the following section we will discuss in greater detail the properties of the rolling 
moment and its implications on the structure and stability of agglomerates. 

§ 4. Mechanical properties of the contact under the influence of rolling 

MOMENTS 

4.1. Static considerations 

Equation (25) shows that the force opposing the attempt to roll one particle over 
the other is a torque that is proportional to the angular displacement. However, as the 
displacement increases there is a limit to the increase in that force. When the line 
connecting the two particle centres has moved a certain critical distance over the 
surfaces, the contact area will start to move, and the force will decrease or drop to zero. 
We will call this critical distance the yield displacement £ yield = /?<5a‘ dd . £ yield should be 
of the order of the distance between atoms in the materials, typically 2 A. Therefore, 
the adhesive contact reacts elastically in a small displacement interval. For larger 
displacements, energy is dissipated and the contact area moves irreversibly over the 
surfaces. 

Due to the simple structure of eqn. (25), the calculation of several mechanical 
properties is quite straightforward. We will illustrate the results in some figures, which 
show the involved forces and other quantities as functions of the surface energy y and 
of the particle radius R (see figs. 3 and 6). 

For the following discussion we will restrict ourselves to the contact between 
identical spheres with radius R\. Then, R = R\/2 and P c = \nyR\. Furthermore, we 
neglect the variations in the rolling moment due to variations in the normal force. These 
enter into eqn. (25) in the factor a m . However, if the normal force stays within the range 
— P c < P < P c , the factor a 1 ' 5 varies only in the range between 0-5-12. Larger normal 
forces are likely to eventually destroy the contact, and are therefore not relevant to our 
problem. Therefore, we assume for most of the rest of the paper that d= 1. 

Since we have neglected the normal load, the surface energy y is the sole material 
property that enters into the equations. In the diagrams, we will therefore usually plot 
y on the abscissa. Every material or combination of materials may then be associated 
with a position on this axis. We have plotted the positions of some materials on the top 
of the diagrams. The material properties used are summarized in the table. The ordinate 
of the plots is labelled by the reduced particle radius. Therefore, for a given system, 
the results can be read immediately from the figures. 

The first quantity of interest is the tangential force that can be applied at the centre 
of the particles in order to exceed the yield displacement. In other words, to which 
limiting force does the contact still react elastically? A force £? yield , applied at the particle 
centres, corresponds to an applied moment at the contact point of Af yield = Q y,eld /?i, thus 
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Plots of some results for the adhesive contact of identical spheres in the y-R\ plane, (a) The 
pull-off force as calculated from JKRS (P c - 3 nyR) (the labels marking the curves denote 
log Pc [N]), ( b ) the yield force required to initiate rolling (see text) normalized by the 
pull-off force Qy^JPc, (c) the minimum radius of curvature ^ min [m] that can be applied 
to a linear chain without leaving the elastic domain, and (d) the number of particles 

that can be in a linear chain of particles in the earth’s gravitational field without leaving 
the elastic domain. Calculated for particles with a mass density of 2 gem -3 . 


0 yietd = 6Ti}’£ yield , 


(26) 


or 



(27) 


The absolute force £) yie,d needed to leave the elastic domain is independent of the 
particle radius, while the force in units of the pull-off force is independent of the surface 
energy y. This behaviour is also reflected in the map of the normalized yield force in 
fig. 3 ( b ). There it can be seen that for micron-sized or larger particles, the force required 
to initiate rolling is much smaller than the pull-off force. Therefore, it should be quite 
simple to roll particles of these sizes over each other, without immediately destroying 
the contact. On the other hand, for very small grains, the force needed to initiate rolling 
becomes comparable to the pull-off force. For particle radii of 1 nm, the forces would 
even be the same — however, we have to bear in mind that at these sizes 
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Material properties. 


Material 

yt (niJ m 2 ) 

E (N m -2 ) 

V 

p (gem 3 ) 

Reference 

Polystyrene 

12 

34(9) 

0*5 

1*04 

[2,3] 

Quartz 

25$ 

54(10) 

0-17 

2*6 

[1,2] 

Graphite 

75 

10(10) 

0*32 

2*2 

[4,5] 

Gold 

790§ 

— 

— 

19*3 

[7] 

Iron 

3000 

2-1(11) 

0*27 

7*7 

[2,6] 


t Surface energy per surface. 

$ Measured for micron-sized particles. 

§ Mean value of the range given in the reference. 

References; [1] Kendall, Alford and Birchall 1987, [2] Anderson 1981, [3] Kendall and 
Padget 1987, [4] Brocklehurst 1977, [5] Zisman 1963, [6] Easterling and Tholen 1972, 
[7] Hodgman 1949. 


the theory of elasticity is no longer applicable. The fact that both forces are comparable 
implies that any attempt to compress an aggregate by rolling the particles over each 
other would be likely to destroy the aggregate. 

Another interesting quantity is the flexibility of particle chains. Many open 
aggregates mainly consist of quite linear chains in a tree-like structure. A way to treat 
this problem is to calculate the minimal radius of curvature that a linear chain of particles 
can adopt, again without leaving the elastic domain. If we bend every contact to its limit 
(which is the angle <5 yield = 2£ yseld //?i), the radius of curvature of the chain is 


$ min 


R] 

^ yield * 


(28) 


a result that is mapped in fig. 3 (c). The minimal radius of curvature is generally large 
compared to the radius of the particles. Even for rather small spheres (lOOnm), the 
minimal radius of curvature is about 50 pm, thus a factor of 500 times larger. This means 
that a chain would have to be a few thousand grains long in order to be elastically bent 
to a circle. Chains shorter than this length would appear as quite rigid structures that 
could not be reversibly bent. Because of this rigidity, fractal ‘fingers’ of aggregates do 
not close up into ring-like structures easily. 

In a laboratory, the aggregate is usually subjected to the earth’s gravity. Another 
interesting question is, therefore, whether an aggregate will collapse under its own 
weight and flatten. We will restrict ourselves to a very simple consideration. We again 
look at a linear chain of N particles bonded by adhesion and fixed at one end to the wall 
(c.f. fig. 4). We make the further simplification that the chain is not bent, but assume 


Fig. 4 



Experiment for the determination of the influence of the earth's gravitational field on the stability 

of an aggregate. 
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that it is a straight line. This is certainly justified by the above discussion of the radius 
of curvature. We may then compute the torque that acts on the key contact of the chain 
(which is the strongest torque anywhere in the chain). It is given by 

= mgRiN 2 , (29) 

where m is the mass of an individual particle, g is the gravitational acceleration and N 
is the number of grains in the chain. Irreversible roiling occurs now when N > A^ grav , 
where 



and p is the specific density of the grain material. We have calculated for 
p = 2gcm -3 and plotted this in fig. 3. The results hardly depend on the density, since 
jVg ra v scales only with the square root of p . It is obvious that for particles larger than 
about 100 pm, even a Single particle is heavy enough to start rolling under the earth’s 
gravity. However, with decreasing particle size, the number of grains that can be 
supported increases rapidly. A chain of approximately one thousand lOOnm spheres 
is still stable under the earth’s gravity. A more careful analysis of the stability of tenuous 
structures (not limited to straight chains of particles) was carried out by Kantor and 
Witten (1984), using a parametrized strength of the contact. The size limits established 
in their study are in good agreement with our results when appropriate estimates for 
the contact strength are used. 

4.2. Dynamic considerations 

4.2.1. Vibrations 

As already noted above, for small angular displacements the moment produced by 
the displacement is proportional to the displacement itself. It should therefore be 
possible to excite vibrations where two particles in contact oscillate around the 
equilibrium state as shown in fig. 5 ( b ). Before we calculate the frequency associated 


Fig. 5 




The two modes of vibration considered in § 5.2. 1 . {a) Vertical vibration: the vertical displacement 
S z is oscillating, when the distance between the sphere centres increases and decreases; 
( b ) rolling vibration: the angular displacement S ay oscillates, when the spheres roll back 
and forth over each other. 
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with these oscillations, we will briefly discuss (as a means of comparison) a different 
type of vibration in which the particles oscillate along the z axis as shown in fig. 5 (a). 

If v z \ and v Z 2 denote the velocities of the particles in the z direction (defined as 
indicated in fig. 5), the equations of motion are 

miV zl =P, (31) 

m 2 v z 2 = -P. (32) 

The change in the vertical displacement is given by 6 Z = v z \ — v z2 , thus 


$ = V Z I “ V Z 2=—P, 
m 


(33) 


where m = ( \fm\ + \lm 2 )~ 1 is the reduced mass of the particles. P is of course not a 
linear function of S z , but we can linearize it in the vicinity of the equilibrium point P — 0 
(see also Tholen (1992)) and find 




6 E*a 0 s 
5~ S 


The frequency of the vibrations is therefore given by 



(34) 


(35) 


If we consider again identical spheres of radius R x and density p and with 
£* = E\!( 1 — v 2 ), we find 


^ 6 . 38 x 108 ( * 

\10 °Nm -2 / VIOOmJm" 2 / V R } ) \ p } K } 


Thus the frequencies are typically of the order of GHz, and almost independent of the 
surface energy of the materials. 

For the angular vibration we neglect all interactions with other degrees of vibration. 
The equations of motion then become 

m\k]cb x =My(d <xy \ (37) 

m 2 k\d >2 = - Myidzy), (38) 


where cu, is the angular velocity of particle i about the y axis (as indicated in fig. (5)) 
and ki is the radius of gyration of the particle ((2/5) 1/2 times the radius for homogeneous 
spheres). The time derivative of the angular deformation is given by cb\ — d) 2 , thus 


CO\ — OJ 2 


J 


1 1 
Im]/: 2 m 2 k 


)My = -jM y 


4 PcR 

I 


S 




(39) 


where I is the reduced moment of momentum of the two spheres. Thus, the system 
oscillates in this degree of freedom with an eigenfrequency of 



( 40 ) 
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Fig. 6 



(a) W 


Eigenfrequencies for the z vibration (a) and the rolling vibration (fc) for two identical spheres 
in contact, plotted as contour lines in a radius-surface tension diagram. The density of 
the materials was assumed to be p = 2 gem _ 3 . In (a) the quantity E* (see text) was set 
to 10 10 Nm -2 . The curves are labelled by log (frequency) in units of Hz. 


Considering two identical homogeneous spheres with radii R\ and a mass density, 
p, we have 

/? = /?,/ 2, *, =k 2 = (H5) m Ru m, = m 2 = (4/3)u/?fc, 



These frequencies as functions of y and R\ are shown in fig. 6. The eigenfrequency 
associated with the rolling motion is usually one to two orders of magnitude smaller 
than the frequency of the vertical vibration. 


4.2.2. Energy dissipation during rolling 

The rolling vibration can only take place if the amplitude is small enough. For 
amplitudes larger than a critical angle <5£ e,d = £ yield //?, the spheres start rolling 
irreversibly over each other. When we assume that after a rolling distance £ the 
contact area readjusts to a symmetric contact with zero torque, the average torque is 
(owing to the linear dependence between torque and £) 

(My) ~ 2P C £ yidd , ( 43 ) 

which is equivalent to an average decelerating force applied at the centre of the rolling 
sphere of 

cyield 


( 44 ) 
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Fig. 7 



Experiment to show the amount of energy lost due to rolling friction. One of the two spheres 
is fixed in an inertial system (hatched area). The other sphere is approaching with a kinetic 
energy W 0 . Upon contact, the energy W 0 is distributed between rotation and linear 
movement. The moving sphere rolls over the fixed sphere until all energy is dissipated. 


In order to visualize this energy loss we will calculate the number of full rotations 
a sphere can roll with a given energy at the start. We consider a sphere colliding with 
a fixed sphere (see fig. 7). We assume that sliding friction is strong enough to entirely 
prevent sliding, so that the surfaces do not move relative to each other. Therefore, the 
initial energy W 0 is distributed into rolling and linear motion: 

Wb = ±mir? + ±/,co?. (45) 


Since no sliding occurs, V\ — coiRi. If we further observe that for a homogeneous sphere 
I — we find that the rotation of the sphere just after contact is 


co(t = 0) = 



(46) 


From then on the equation of motion a)) is given by 


d>i = ” 


My 

Tr 


(47) 


where M y is given by eqn. (43). Integrating eqn. (47) we can find the time when the 
rotation of the sphere stops, and calculate the number of rotations performed. The 
result is 


1 Wo 

rtro11 1471 P c £* icld ' ^ 

Thus, the number of rotations is proportional to the initial energy and inversely 
proportional to the pull-off force. This seems to indicate that rolling is much more likely 
in collisions between grains with a small pull-off force. However, it is also to be 

recognized that grains with a small pull-off force will only stick at small impact 

energies. We should therefore measure the energy W 0 in units of an energy typical for 
adhesion. Such an energy is given by P c 8 0 . It has been shown by Chokshi, Tielens and 
Hollenbach (1993) that this energy is of the order of the maximum energy that can be 
dissipated in a slow collinear collision between two spheres, so that they will stick 
together. Then, eqn. (48) becomes 

_ 3 / y N 273 R m / Wo \ 

” ro11 14(9 jt) 1/3 U*j £ y,eld \P C <V‘ (49) 

Thus the material properties enter only in the form of the ratio E*/y and we may 
plot M r0 H in a map as shown in fig. 8. Even in a collision with relatively large energies 



Resistance to rolling 


797 


-2 


-4 

£ 

DC* 

? -6 


-0 

9 10 11 12 13 

log E*/y [m-‘] 

Plot of the number of rotations a mobile sphere rolls over a fixed sphere with an initial tangential 
kinetic energy P c Sq> The curves are labelled by n r0 n/ ( Wo/Pc^o)- 1° the region to the lower 
right, the rolling distance is smaller than £ y,etd . 

(Wo 2=55 Pc<>o), spheres with a large ratio £*/y do not roll at all, if they are smaller than 
about lOnm. Larger radii, smaller elastic moduli, or larger surface energies are 
favourable to more rotation — however, the rotation is usually less than one complete 
rotation if the radius of the spheres stays below 100 nm. For bigger spheres, it becomes 
quite possible that rolling persists for a full rotation or more. This means that in a 
collision, the contact would not remain at the point of first contact. Rather, the spheres 
will roll over each other. If the fixed sphere is part of an aggregate, the mobile sphere 
will probably roll until at least a second point of contact is established. Therefore, the 
resulting aggregates will be considerably more compact. 

§5. Conclusions 

In this paper we have discussed the resistance of a JKRS contact to rolling. This 
was done by recognizing that, due to the atomic structure of the surfaces on contact, 
the contact area will not be symmetric around the contact point. This asymmetric contact 
area leads to a net torque force opposing the attempt to roll one body over the other. 
We have calculated the torque of the asymmetric contact area by an approximation in 
which the area was divided into two half cycles of different size. This approximation 
was justified by an exact solution for the two-dimensional case. However, it would be 
very interesting to see an exact analytical solution for the three-dimensional case. 

Since the maximum asymmetry of the contact area is given by the typical length 
in the atomic grids of the bodies irreversible rolling starts only when the bodies roll at 
least this critical distance. For smaller distances, the contact reacts elastically with a 
force proportional to the displacement and to the pull-off force. As long as the contact 
is not subjected to normal forces, the resistance to rolling is dependent on the surface 
energy alone. Normal forces introduce a very weak dependence upon Young's modulus. 

We have also discussed some implications of our results, many of which should be 
testable by experiments. The results indicate a critical bending angle below which the 


Fig. 8 
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bending will be elastic. The critical angle is quite small, so that long linear chains of 
particles would appear rather rigid. We also showed that in the case of very small 
particles the force needed to bend a chain can be comparable with the pull-off force. 
Therefore, aggregates of very small grains would be very stiff. Forces applied to these 
aggregates are likely to break them instead of deforming them considerably. 

Finally we have shown that it should in principle be possible to excite vibrations 
not only in the vertical degree of freedom, but also in the rolling degree of freedom. 
This should be visible in vibrations of fractal aggregates as well. 
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APPENDIX A 

Comparison of exact and simplified calculation of the rolling moment of an 

ASYMMETRIC ADHESION PROFILE 

In this appendix we will calculate the pressure distribution and the associated rolling 
moments of an elastic cylinder in contact with a rigid plane. In this 2D case, it is possible 
to calculate the pressure distribution in an asymmetric contact analytically. This will 
enable us to compare the rolling moment obtained from a simplifying assumption as 
used in §4.1 with an exact solution. 

Although we here describe only a contact between an elastic cylinder and a rigid 
plane, the results can be immediately generalized to the case of two elastic aligned 
cylinders by just replacing the cylinder radius R and the elastic constant E* = £7(1 — v 2 ) 
with reduced values. As usual, E is the elastic modulus of the material and v is the 
Poisson’s ratio. 

Both the symmetric and the asymmetric contact are shown in fig. 9. Our treatment 
of this 2D contact problem follows the description by Johnson (1989). 


A. 1 . The symmetric pressure distribution 
When the cylinder is in contact with the plane, the deformation u z of the surface 
of the cylinder is given by 


u z (x) 



— a < x < a. 


(50) 


Fig. 9 



2a 


Symmetric and asymmetric contact of a cylinder with a rigid plane. 
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where we put the origin of the coordinate to the same x coordinate as the centre of the 
cylinder. In the case of specified deformation in the loaded region, the pressure 
distribution p{x,a) can be calculated from the following integral equation (see, for 
example, Johnson (1989)): 


gCMjO 

X~S 


d.v = — 


nE du z 
2 dx 


(51) 


This is an integral equation of the type 

f =*(,>. 

J~ a x-s 


(52) 


The general solution of this equation can be obtained by applying the finite Hilbert 
transformation which leads to 


P(x, a) 


1 


JtV -* 2 ) 1 


[f 


(a 2 - 5 2 ) ,/2 g(s) 


ds + C 


x — s 


]■ 


(53) 


where C is some constant of integration. In the following we add a subscript v to the 
pressure distribution in order to indicate that this is the symmetric case. The application 
of eqn. (53) to eqn. (51) yields 


Ps(x , a) = 


1 



itE sfa 2 -* 2 ) 1 ' 2 
2 R x — s 


E 

2/?(a 2 — Jt 2 ) 1 ' 2 



( a 2 -s 2 ) m 


x — s 


ds j. 

(54) 

dsj. 

(55) 


The integrand is singular at x = s and we have to take the Cauchy principal value, which 
is given in eqn. (81). Inserting this result and rearranging yields 




0 a 2 — x 2 ) 1 ' 2 + (C — a 2 ! 2) 


(a 2 — x 2 ) l/2 J 


(56) 


The first part of the pressure distribution is the usual Hertz solution for two cylinders 
pressed into contact. In this part of the solution the pressure falls to zero as x goes to 
a. The second part which is proportional to l/(a 2 -x 2 ) u2 is the one-dimensional 
equivalent to the solution of Boussineq’s problem which shifts the contact area 
vertically by a certain amount and forms the neck which is typical for adhesive contacts. 

The total force per unit length P is: 

P=f p,(x,a)dx (57) 


With the integrals given in eqns. (84) and (83) this becomes 


P = 



(59) 


The rolling moment M y = Jjrp s (jt,a)dx: of this pressure distribution is of course zero. 
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A. 2. The asymmetric pressure distribution 
We now consider a contact where the contact region is of the same size as in the 
symmetric case, but is shifted with respect to the centre of the cylinder by an amount 
£, (see fig. 9). Then, the deformation again follows eqn. (50) but the contact strip is a 
shifted interval: 


Uzto = 5z~ 


X>_ 

2 R' 


— a — £ < x < a — 


(60) 


For the calculation of the pressure distribution it is convenient to have the contact 
area symmetric about the origin. In order to do so we introduce the variable 

X = x+£. (61) 

Then the deformation can be written as 

u z (X) = d z - (<X ^ , -a<X<a. (62) 

2K 


Introducing this result into eqn. (51) yields for the asymmetric pressure distribution 
P*(X, a) 


P*(X,a) = 


2 R(a 2 - X 2 ) 


2x1/2 


C + - I 


If 


X-s 


r 

n J 


(a ~ s 2 ) 1 
X-s 


ds . 


(63) 


With eqns. (81) and (82) this can be written as 

•> ■ - 5 [*■ ~ x ’ ya + ( c - T + ») ra*]- 

Thus the pressure distribution differs from the symmetric one only by an additional term 
which is proportional to X and, therefore, not symmetric. As we will see, this term is 
responsible for the moment associated with this pressure distribution. 

The total force associated with this distribution is again (see eqn. (57)) 


P = 



(65) 


and the rolling moment about the y axis which intersects the x axis at X = 0 is 
given by 


M 


x = 


0 



Xp a (X,a)dX 


( 66 ) 


e r K x 2 
2RJ-J (a 2 — X 2 ) 1/2<iX 


(67) 


7r a 2 E 


t. 


4 R 


( 68 ) 
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We still have to transform this result back in the original coordinate system. The moment 
about the y axis that intersects the x axis at x = 0 can be calculated by 

M* = ( = mX = °-ZP, (69) 

and thus 


M r 0 ~~^ 2 + 2O. (70) 

This is the exact result which we will compare with the approximate calculation 
performed in the following section. 


A. 3. Simplified calculation of the rolling moment 
The moment of the asymmetric distribution can also be approximately calculated 
in the following way. We assume that the x< 0 half of the contact area has a pressure 
distribution equal to the symmetric distribution p s (x, a 4- f) with contact strip half-width 
a + £, while the x > 0 part of the contact area has a distribution equal to a symmetric 
distribution p s (x y a - £)» thus with respect to a contact half-width of a - The overall 
distribution calculated in this way is therefore not continuous at x = 0. 

The moment associated with this distribution is 


«-r.. 


xp s (x 9 a + £) dx + 


- £ 

xp s (x, a 

0 


-$)dx 


r « + £ r<* ~ c 

— I x/? s (x, a -f f ) dx + x/> s (x, a — £) dx 

Jo Jo 


= - , + M {U£} = 

lYt y,a + £ ' 1 Vi y, a — ^ 


0 / 2 ) 


(71) 

(72) 

(73) 


Thus we have to solve integrals of the kind 


W v '3= f xp s (x,. 
Jo 


s/)dx 


E 
2 R 




x(stf — x^ 

The required integrals are given in eqns. (85) and (86). The result is 




(74) 

(75) 

(76) 


The total moment is now 


M y = - M 


(1/2) 

y. a + l, 


+ m {U2) . 

1 lvl \\a - 


which, to first order in £fa, is 


My 


1 E 

2 R 


i(a 2 - 2 C). 


(77) 


(78) 


Comparing this result to eqn. (70), we can see that eqn. (78) is smaller by a factor 
7t/2, but otherwise reproduces correctly the dependencies upon the properties of the 
cylinder and upon the displacement c;. This seems to be a very satisfactory result. It 
justifies the use of the approximate calculation of the rolling moment. 
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APPENDIX B 

Integrals used 

The following integrals have been used in the paper. We write them in terms of two 
abbreviations: 

W ] (r,a) = (a 2 -x 2 ) u \ (79) 

W 2 (r , , a) = ( 1 - r z /a 2 )' n . (80) 

B.l. Cauchy principal values of some singular integrals important in the 
two-dimensional problem 


f a s(a 2 — s 2 ) 172 2 2 

P. V. — — ds = Kix 1 - \a 2 ), 

J-a X-S 

(81) 

r (a 2 — s 2 ) m 

P.V. — d.v = jut. 

J -a X - s 

(82) 

B.2. Non-singular integrals 


r i . 

dx = 71, 

J^Wfr^a) 

(83) 

f W\(r, a)dx — \na 2 , 

J — a 

(84) 

f xW\(r,a)dx= £a 3 , 
Jo 

(85) 

f * d ,-a. 

Jo W,(r,a) 

(86) 

ra ^ 

SW 2 (r,a)dr = —a\ 
Jo lb 

(87) 

f E . n , 

d r = —ar. 

Jo Wfr.a) 4 

(88) 
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